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ABSTRACT 

We present self-consistent star formation rates derived through pan-spectral analysis 
of galaxies drawn from the Galaxy and Mass Assembly (GAMA) survey. We determine 
the most appropriate form of dust obscuration correction via application of a range of 
extinction laws drawn from the literature as applied to Ha, [On] and UV luminosities. 
These corrections are applied to a sample of 31 508 galaxies from the GAMA survey 
at z < 0.35. We consider several different obscuration curves, including those of Milky 
Way, Calzetti (2001) and Fischera and Dopita (2005) curves and their effects on the 
observed luminosities. At the core of this technique is the observed Balmer decrement, 
and we provide a prescription to apply optimal obscuration corrections using the 
Balmer decrement. We carry out an analysis of the star formation history (SFH) 
using stellar population synthesis tools to investigate the evolutionary history of our 
sample of galaxies as well as to understand the effects of variation in the Initial Mass 
Function (IM F) and the effects this has on the evolutionary history of galaxies. We 



find that the Fischera & Dopita (20051 obscuration curve with an R v value of 4.5 
gives the best agreement between the different SFR indicators. The 2200 A feature 
needed to be removed from this curve to obtain complete consistency between all SFR 
indicators suggesting that this feature may not be common in the average integrated 
attenuation of galaxy emission. We also find that the UV dust obscuration is strongly 
dependent on the SFR. 

Key words: galaxies: evolution - galaxies: formation - galaxies: general - galaxies: 
initial mass function 
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1 INTRODUCTION 

The absorption of light by dust in galaxies presents a major 
challenge in measuring accurate star formation rates (SFRs) 
in galaxies. Understanding star formation is essential to our 
understanding of the evolutionary history of the Universe 
on both local and global scales. SFRs are commonly derived 
from Ha and ultraviolet (UV) luminosities, both of which 
are significantly affected by dust obscuration. The basic ap- 
proach to correcting for this relies on having an observed 
constraint on the level of the attenuation, such as the Balmer 
decrement ( jSeaton|1979l|Cardelli et al.|1989| |Calzetti|2001 1 



or the /3-parameter (Calzetti, Kinney 



Storchi-Bergmann 
, |1999| |Adelberger ~ 
Steidel 2000 I together with a model for the wavelength- 



1994; IMeurer, Heckman and Calzettil 



Cardelli et al. 



dependence of the at tenuation (|Seaton||1979 
1989) |Calzetti|[200l| [Pierini et al.||2004[ |Tuffs et al.||2004| 
Fischera fc Dopita||2005| ~ 



Dust has the dual effect of not only reducing the overall 
luminosity observed from a galaxy but also imprinting wave- 
length dependant spectral variations ( |Calzetti|1997a| ). Dust 
absorbs short-wavelength radiation and re-emits this energy 
at longer wavelengths. Dust also plays a significant role in 
the fluid dynamics, chemistry and star formation of galax- 



ies (Draine et al. 20071. Even though dust plays a crucial 
role in galaxies, quantifying its effects and physical proper- 
ties remains challenging. Understanding the roles played by 
inclination angles, optical depths and galaxy morphologies 
when coupled with dust has been even more difficult, and 
only recently have these effects been quantified by |Tuffs et| 
aLlpOOil). 



Correcting for the effects of dust attenuation is usu- 



ally done using semi-empirical models (see Cardelli et al 



1989 Byun et al. 1994 Calzetti 20011 but there are ra- 



diative transfer (RT) models developed for attenuation cor- 



rections (Xilouris et al.||1997 1998 1999 Witt 1992 Tuffs 



et al. 20041. Even though modelling techniques can give 



more accurate attenuation corrections than semi-empirical 
models, they require more observational information about 
the galaxies in order to produce accurate corrections on 
an object-by-object basis. However, the RT techniques have 
proven extremely powerful when utilised to create libraries 
of model simulations ( Tuffs et al.|2004 \ for use in statistical 
analysis of large statistical samples ( jDriver et al.| :2009). 

An important method for providing accurate SFRs is 
the spectral energy distribution (SED) modelling technique. 
There are two main types of SED techniques: template fit- 
ting SEDs ( [Calzetti et al.||2000| |Salim et al.||2007 | 



SED techniques ( Silva et al.||1998 Popescu et al. 



and RT 
20001. In 



particular the accuracy of the RT SED techniques is high, 
but it tends to be complex and it requires knowledge of the 
SED over the whole electromagnetic spectrum, from UV to 
far-infrared (FIR) and sub-millimeter wavelengths, which is 
not always available. Due to this complexity the use of RT 
techniques are less favoured than empirical methods for dust 
correction. 

Balmer decrements are commonly used as an indicator 
of obscuration corrections and are usually applied in con- 
junction with an extinction law such as those for the Milky 
Way (MW) ( |Nandy et al.|1975| |Seaton|1979| jCardelli et al 



1989|), Smal l Magellanic Cloud ( |Prevot et al.|1984||Bouchet 
et al.|| 1985| SMC), Large Magellanic cloud (|Fitzpatrick et 



al.|1986[ LMC), th at observed for M31 (jBianchi et al.|1996[ ), 



starburst galaxies ( Calzetti|19 97a, 20011 and other galaxies 



(Fischera & Dopita 2005 FD05 hereafter). Any extinction 
law can be used to correct either nebular emission lines or 
continuum emission. As a consequence of how these curves 
are derived, though, some curves are preferred over others to 
correct for each type of emission, as the processes that result 
in these two types of emission are different and it has been 
shown that nebular lines are usually obscured more than 



the continuum emission l 


Fanelli et al. 


1988 


Calzetti, Kin- 


ney & Storchi-Bergmann 1994 , Calzetti 


1997 


)||Mas- Hesse & 



Kunth 19991. A physical interpretation for this is provided 



by Keel ( 1993| ) and Calzetti, Kinney & Storchi-Bergmann 



( 1994 1 who argue that the ionizing hot (young) stars which 
produce the nebular lines are found in or close to the dusty 
molecular clouds from which they were born, while the UV 
continuum is a product of older stars that have over time 
moved away from the dust clouds in which they formed, or 
have destroyed the dust in situ revealing the cluster of stars. 

The obscuration curves are also dependent on their en- 
vironments, in particular the metallicity of the host galaxy 



( jCardelli et al.|1989[ ). For instance, the LMC, SMC and MW 
curves have a large variation in values at the FUV end of the 
spectrum due to the different metallicities of the three galax- 



ies and grain size distribution (Calzetti, Kinney & Storchi- 
|Ber gmann|1994| |Ca lzctti 1997b) in the interstellar medium. 
These LMC, SMC and MW obscuration curves only take 
into account the dust between the observer and the star but 
in the case of external galaxies the dust geometry is more 
complicated. Not only does the dust between the observer 
and the stars of the galaxy need to be taken into account, 
but back-scattering of light into the line of sight from dust 
in other regions of the galaxy must also be considered and 
the dust is never uniformly distributed. The Calzetti extinc- 
tion curve attempts to consider all these issues by folding a 
variety of obscuration effects into a single expression. 



Meurer, Heckman and Calzetti ( 1999 1 show a relation- 



ship between the ratio of far-IR (FIR) and UV fluxes and 
the UV spectral slope /3 for a sample of starburst galaxies. 
Because this technique relates the FIR and UV radiation 
emitted from galaxies it can be a powerful tool in recover- 
ing the UV radiation lost due to the dust, regardless of the 
geometry of the dust. The /3 parameter is based on the rela- 
tion between the gradient in the UV obscuration curve and 
the UV wavelengths ( Meurer, Heckman and Calzetti||1999 



Calzetti, Kinney fc Storchi-Bergmann 1994; Adelber ger fc| 



Steidel 



2000) 



However, this method is only recommended 
for starburst galaxies as quiescent galaxies tend to deviate 
from the total FIR to UV luminosity ratio and UV spec- 
tral slope relation seen for starbursts (Kong 20041. The ac- 
curacy of this method is highly dependent on the method 
employed to measure the UV spectral slope. |Calzetti, Kin-| 
ney fc Storchi-Bergmairn] ( |1994 1 used 10 bands along the 



observed frame UV continuum spectrum avoiding all large 
scale features that deviate from the trend of the slope. With 
large scale surveys such as GALEX this level of accuracy is 
unattainable, leading to less accurate measurements of the 
slope and potentially flawed dust corrections. 

Once suitable dust corrections are in hand we can de- 
rive accurate SFRs and examine the relationship between 
galaxy properties and the dust content in a galaxy. SFRs 
are related to the luminosity through a linear scale factor, 



Obscuration in galaxies and SFHs 3 



determined with the assumption of a constant stellar initial 
mass function (IMF). These conversion factors will change 
if the IMF is varied. One recently suggested possibility is an 
evolving IMF ( |Wilkins et al.|2008a|b| Gunawardhana ct al. 



m prep 



which would lead to an evolving SFR conversion 



factor. Pnarnrn-Altenburg, Weidner & Kroupa (20091 pro- 



posed an Integrated Galaxy IMF (IGIMF) that combines 
the IMFs of all young-star clusters to form a galaxy-wide 
IMF. This was developed to account for the inconsistencies 
in current IMFs which are based on isolated stellar clusters 
and then applied on galaxy wide scales. One of the incon- 
sistencies that the IGIMF accounts for is the discrepancy 



between FUV and Ha derived SFRs at low SFRs (Sullivan 
et al.|2000||Lee et al.||2009 1 



We compare and contrast SFRs corrected for dust using 
different obscuration correction methods, identifying an op- 
timum approach, and then use that in a preliminary investi- 
gation of SFR histories of galaxies. The SFRs are compared 
with theoretical evolutionary paths to better understand the 
star formation history of the galaxies in our sample. 

In §2 we present details of the data used. In §3 we 
present a prescription for the derivation and the application 
of the obscuration corrections. §4 compares Ha, FUV and 
[On] derived SFRs. §5 examines the evolution of our SFRs 
by comparing these SFRs to evolutionary synthesis models. 
We also compare the evolutionary paths of the SFRs with 
the predictions of the IGIMF theory. §6 is the summary 
and conclusion of our findings. The assumed cosmological 
parameters are: 

Ho =70kms" 1 Mpc" 1 , Q M = 0.3 and D. A = 0.7. All 
magnitudes are in the AB system. 



2 DATA 



We use data from the Galaxy and Mass Assembly (GAMA) 
survey ( jDriver et al.|2009| ). GAMA is a multi-band imaging 
and spectroscopic survey covering «144 square degrees of 
sky in three 12° x 4° regions ( Robotham et al.|2009 Baldry 



et al. 2009). The spectroscopy comes from the AAOmega 



spectrograph ( Sharp et al. 2006 1 at the Anglo Australian 



Telescope (AAT). 

Th e UV data was obtain ed from the GALEX-GAMA 
survey (Seibert et al. in prep. I. A signal-to- noise ratio cut of 



2.5 was applied to the UV data. There are 111521 GALEX 
sources within the GAMA regions of which 110443 matching 
sources were found using a matching radius of 6". Such a 
large matching radius is is necessary due to the increased 
positional uncertainty inherent in the larger GALEX PSF. 
SDSS and spectral line data were matched to GAMA sources 
using a matching radius of 2.5" resulting in a final sample of 
47269 objects. After removing AGN (using the prescription 
of Kewley, 2006) and sources without emission lines, r-band 
magnitude and redshift data the final sample size contained 
31508 galaxies. All remaining galaxies contain both Ha and 
H/3 fluxes. The final sample has a GALEX FUV magnitude 
range of 14.7 < tufuv < 25.3 and a redshift limit of z < 
0.350 due to the requirement for Ha being in the observed 
spectral range. 



3 DUST OBSCURATION CORRECTIONS 

Correcting galaxy emission for dust requires a case spe- 
cific approach, that is nebular and continuum emission from 
galaxies require two different treatments. At the core of these 
obscuration corrections are the obscuration curves and the 
Balmer decrement. We apply several different obscuration 
curves to examine their effects on continuum (UV) and neb- 
ular (Ha and [On]) emission. 

3.1 Continuum emission corrections 

Applying obscuration corrections for continuum emission re- 
quires two parts, the reddening suffered by the stellar con- 
tinuum (E(B — V) CO nt) and an applicable obscuration curve 
(fc(A)). We use the standard form 



L o 10 



0AE(B-V) aont k(X) 



(1) 



to apply our dust corrections. Li and L a are the intrinsic and 
observed galaxy luminosities respectively. To obtain k(\) 



for the UV stellar continuum, the Calzetti ( 2001 1 obscura- 



tion curve and the FD05 obscuration curves can be applied. 
When using the [Calzetti (2001) curve fc(A) must be divided 
by 0.44. The FD05 curves can also be applied to correct for 
the obscuration of nebular emission. The UV luminosities 
were derived using the UV fluxes from the GALEX survey 
and the redshifts obtained from the GAMA survey. 

E(B — V) con t is derived from the reddening in the ion- 
ized gas, E(B - V) ga s- E(B - V) cont = O.UE(B - V) ga s 
( Calzetti|200lj ) and 



lot 



E(B - V) gas = 



(fo/ 2 - 86 ) 
0A[k{H/3) - k(Ha)] 



(2) 



where k(H(3) and k(Ha) are the obscurations of the nebular 
emission at the Ha and H/3 wavelengths derived from a MW 



obscuration curve ( Seaton 1979 Cardelli et al. 1989 \ or a 



theoretically modelled curve (FD05) as in this paper, fna 
and /h/3 are the stellar absorption corrected but not dust 
corrected Ha and H/3 fluxes. 

fna and fs/3 were corrected for stellar absorption using 



the formalism outlined in Hopkins ct al. ( 2003 ) stated below 



EW + EW C 



EW 



fo 



(3) 



where f a and f s are the observed and stellar absorption 
corrected fluxes respectively. EW is the equivalent width of 
the line being corrected and EW C is the correction for stellar 
absorption taken to be 0.7 A which has been shown to be a 
reasonable assumption ( jGuna wardha na" et al. in prep.[ ). 

Assuming a Case-B recombination with a density of 
100 cm -3 and a temperature of 10 000 K, the predicted ratio 
of f Ha to f Hp is 2.86 ( |Osterbrock[T989| . All Balmer decre- 
ments below 2.86 were set equal to 2.86 as suggested by |Kew-| 



ley ( 2006 1. We applied k-corrections to the observed GALEX 
UV magnitudes, using KCORRECT.V4.1.4 ( [Blanton et al.| 
2003 1 to infer the rest frame magnitude in each GALEX 
band, with effective wavelengths of 1528 A and 2271 A , be- 
fore the dust correction process. 

It is important to note that correcting continuum emis- 
sion for dust using this technique requires both MW and 
continuum obscuration curves. MW obscuration curves such 



as Seaton (1979); Cardelli et al. (19891 are required for the 
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derivation of E(B — V) CO nt while a continuum obscuration 



curve such as Calzetti (20011 provides k(X) 



A theoretically modelled curve such as FD05 is appli- 
cable to both the nebular and continuum components. The 
reason for this is that unlike the other extinction correc- 
tion curves, these curves were derived from dust models 
that do not require assumptions about the effects of the 
emission sources of the galaxy. The model is based on the 
inferred physical properties of the turbulent density struc- 
ture (FD05). The FD05 curves are described for a range of 
R v values and we examine the effects of different R v values. 
R v — Ay /E(B — V) ga s, where Av is the extinction suffered 
in the rest frame V-band. 



3.2 Nebular emission-line corrections 

To correct Ha emission for dust, the Ha luminosity needs 
to be calculated first. Ha luminosities corrected for stellar 
absorption, taking into account aperture corrections, can 
be obtained using the formalism outlined in |Hopkins et al7| 
(2003], 



L Ha = (EW Ha +EW c )W 



-0.4(M r -34.10) 



3 x 10 is 
[6564.61(1 + z) 



•(4) 



M r is the fc-corrected absolute r-band AB magnitude. The 
last factor converts units from W Hz -1 to WA" 1 and 
6564.61 A is the vacuum wavelength of Ha. 

Aperture corrections account for the fact that only 
a limited amount of emission from a galaxy is detected 
through the 2" diameter AAOmega fiber. The aperture 
correction implicitly assumes that the emission measured 
through the fiber is characteristic of the whole galaxy and 
that the star formation is traced by the r-band continuum 
emission | |Hopkins et al.|2003| ). 

Obscuration corrections are then applied to the Ha lu- 
minosity as follows, 



Li = Li 



fHct/ f. 

2.86 



116 



(3) 



The exponent in the above equation is derived from the 
FD05 (R v — 4.5) extinction curve and will vary according 
to the (MW) obscuration curve. This form of obscuration 
correction has the same end result as 



U = L Ha lQ 0AE( - B - V) ^ 



fc(A) 



(6) 



In this case k(X) can be derived from a MW obscuration 
curve such as |Seaton| ( [19791 ); [Cardelli et aL l ( fl989j ) or a the- 
oretically modelled curve such as FD05. 

An analogous method can be used to derive [On] lumi- 
nosities but the stellar absorption correction (EW C ) is not 
required. M r is replaced by M u and the Ha wavelength is 
replaced by the [On] doublet effective wavelength, 3728.30 A 
( Hopkins et al.|2 003). As in the case of Ha, aperture correc- 
tions are included but obscuration corrections are not. 



'[On] 



-0.4(M„-34.10) 



3 x 10 18 



[3728.30(1 + 2) 



(7) 



[On] emission can be used as an SFR estimator, par- 
ticularly at high redshifts ( jGilbank et al.||2010[ ). The [On] 
emission must be corrected for obscuration based on the 



obscuration at the wavelength of Ha so that the SFR cali- 
bration derived for Ha luminosities (calibrated to [On]) can 
be used ( |Hopkins et al.||2003| |Kennicutt|1992l |1998| ). 



3.3 (3 Parameter 



The f3 parameter (Meurer, Heckman and Calzetti 1999 



Adelberger & Steidel 2000 1 has been proposed in order to 



provide obscuration corrections for starburst galaxies when 
the Balmer decrement is not available. It is derived from the 
relationship between the UV spectral slope and UV to FIR 
flux ratio. 

The f3 parameter is the UV spectral slope determined 
from a power-law fit to the UV continuum of the form: 



h oc A f 



(8) 



where f\ is the flux density per wavelength interval and A is 



the central rest wavelength ( Meurer, Heckman and Calzetti 
|1999[ ) . The relationship between the j3 parameter and the ra- 
tio of far-infrared and UV fluxes is used to account for dust 



in galaxies. Calzetti, Kinney & Storchi-Bergmann ( 1994 1 



recommend observing the UV continuum bands that avoid 
irregularities along the UV spectral slope that would alter 
the gradient of the slope that is measured. In our case only 
2 UV bands are available, those of far UV (FUV) and near 
UV (NUV) from the GALEX satellite that span 1400 A 
- 1800 A and 1800 A - 2800 A with effective wavelengths 
Apuv = 1528 A and 2271 A respectively. 

Due to the poor sampling of the underlying UV spec- 
trum by the broad GALEX UV bands, and significant con- 
tamination by local structure across these bands the accu- 
racy with which the UV slope can be determined is compro- 
mised. Kong et al. (2004) state that this definition of /3 does 
not apply to galaxies for which only multi-band ultraviolet 
imaging is available, as in the case with GALEX MIS Survey. 
Instead Kong et al. (2004) recommend the definition, 



P = 



logft 



log fa 



logXj 



logXp 



(9) 



where Xfuv = 1528 and Xnuv = 2277 are the effective 
wavelengths of the far and near UV filters of the GALEX 
satellite. }fuv and }nuv are the mean flux densities per 
unit wavelength through these filters. We aim to compare 
the the dust corrections from this technique with corrections 
made using the Balmer decrement. This would allow us to 
observe the effect that the ft parameter will have on dust 
corrections in our sample of largely quiescent galaxies. 

The UV luminosities are corrected for starburst galaxies 
using the relation ^4i6oo = 4.43 + 1.99/3 from Meurer, Heck-| 



man and Calzetti (1999|l, where ^4i6oo is the attenuation at 
1600 A: 



U = 10' 



0.4(4.43+1.99^) 



L . 



(10) 



where Li and L are the FUV dust corrected and uncor- 
rected luminosities. 

We apply this formalism to our sample of largely qui- 
escent galaxies to test the extent to which f3 and hence the 
corrected luminosities and SFRs are overestimated as Kong 
et al. (2004) have suggested. 
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4 COMPARISON OF STAR FORMATION 
RATES 

Star formation rates can be calculated from scaling factors 
applied to the luminosities. The scaling factors are depen- 
dent on the IMF used (|Salpeter|1955| |Kennicutt|1983[ ). We 



use a Baldry & Glazebrook ( 2003 1 IMF with a mass range 
of O.IMq to 120M Q . The effects of apply ing other IMFs 
| Salpeter|1955 Kennicutt 1983; Scalo 1998) are explored in 
§ 4. For Ha, FUV and NUV luminosities the SFR conversion 
factors were calculated using a population synthesis model 



( Fioc fc Rocca-Volmerange||1997| PEGASE.2). 

The [On] SFR conversion factor was derived from the 



flux ratio of [On] and Ha, where -Fjon] I^Ba = 0.23 (Hop- 
kins et a l. 1 2003 1 and the Ha SFR conversion factor ( Kenni- 
cutt |1992|. A stellar absorption correction of 0.7 A 



( Gunawardhana et al. in prep. I is applied before this ratio 



is obtained. The luminosity to SFR calibration for [On] can 
be calculated in this way because the SFRs come from us- 
ing the [On] emission as a proxy for the Ha emission. SFRs 
derived from [On] luminosities are based on the fact that 
there is a good correlation between observed [On] and Ha 
line fluxes ( |Kennicutt|T992} [1998] |Hopkins et al.|2003[ ). 

For FUV, NUV, Ha and [On] luminosities, SFRs can be 
calculated using linear scale factors, derived using PEGASE, 
as follows, 



SFRpwiMp, yr" 1 ) = Lfw^_ 

V (jy ' 1.64 x 10 21 W Hz" 1 



SFRnuv(Mq yr" 1 ) 

SFR [On] (A/ yr- 1 ) = 
SFRhJMq^H 



1.56 x 10 21 W Hz" 

L [0h] 
~ 7.967 x 10 33 W 

Lh<x 

3.464 x 10 34 W 



(11) 
(12) 

(13) 
(14) 



These SFR conversion factors are consistent with those 



of |Kennicutt|(|1998|> for Ha, |Hopkins et al.| ( |2003[ ) for [On] 
and |Madau et aL ( 1998 1 for UV, when taking into account 
the choice of IMF ( jBaldry fc Glazebrook(20(53| . 

Figure [l]a shows the FUV SFRs as a function of Ha 
SFRs with both indicators uncorrected for dust obscura- 
tion, and Figure [l]b shows the same but with only the Ha 
SFRs corrected for dust obscuration. This shows the effect 
that obscuration corrections have on the luminosity and how 
the observed and intrinsic luminosities compare. In Figure [l] 
when dust corrections are applied (in this case only to the 
Ha), the SFR values of high SFR objects increase by ~2 
orders of magnitude. At low SFRs the distribution remains 
close to the one-to-one line. From similar analyses of uncor- 



rected FUV and Ha SFRs by Sullivan et al. (2000) and Lee 



et al. (20091, it is clear that our data follow similar trends. 



Even at high SFR, it is clear that our SFRs are a continua- 



tion of the same trends observed in Sullivan et al. ( 2000 1 at 
lower SFRs. 

Figures[2](FUV) and|3](NUV) show the relationship be- 
tween Ha, FUV and NUV SFRs after the luminosities have 
been corrected for dust obscuration by various obscuration 
curves. Figure[4]shows a similar comparison between Ha and 
[On] SFRs. The obscuration curves are compared in Figure 
5. There is an overestimation of the UV SFRs when only 



UV corrections (Figures [5]a and [3]a). The overestimation 
is larger at higher SFRs. A similar overestimation is seen 
when only the Calzetti| ([2001 ) curve is used to correct both 
Ha and UV SFRs (Figures |2jb and[3]b). The overestimation 
is constant for all SFRs. 

Figures[2]c and[3]c show SFRs after the dust corrections 
were carried out on the appropriate type of emission. This 
means that when dust correcting FUV and NUV luminosi- 
ties a MW obscuration curve ( |Cardelli et al.|1989[ ) was used 
for the nebular emission parts of UV dust correction and a 



Calzetti (20011 curve was used for the continuum parts of 
UV dust correction as outlined in §3.1. The Ha luminosi- 
ties were corrected using the same MW obscuration curve 
( Cardelli et al.|1989[ ) that was used for the nebular emission 
correction part of the UV dust corrections. There is still an 
overestimation in the UV and this overestimation is larger 
than that in panels (a) and (b) in Figure [5] and panel (b) in 
Figure [3] The overestimation is larger at the high SFR end 
for both FUV and NUV. 

In panels (d), (e) and (f) in Figure [2] the same FD05 
curve (i.e. the same R v value) is used for nebular and con- 
tinuum parts of UV dust correction and for the dust correc- 
tion of Ha for each panel. There is good agreement between 
Ha and FUV derived SFRs for an R v value of 4.5 (Fig- 
ure [5]e). This value also agrees with the R v value suggested 
by Calzetti (2001) as well as comparisons made between be- 



tween FD05 curves and the |Calzetti| ( |2001[ | curve (Fischera, 
|Dopita fc Suthe rland 2003]). For R v values less than 4.5, the 
FUV SFRs are overestimated (Figure [2]d, R v — 3.5) and 
for values above 4.5 the FUV SFRs become underestimated 
(Figure [2]f, R v = 5.5), with respect to the Ha SFRs. Both 
Ha and FUV luminosities were corrected using the same 
FD05 curve. 

The corrections applied to the NUV SFR using the 
FD05 extinction curves only lead to an agreement with Ha 
SFRs at an R v value of 6 (Figure[3]f). Applying an R v value 
of 4.5 that is consistent with FUV dust corrections overesti- 
mates the NUV SFRs, as shown in Figure [3]d. R v values 
higher than 6 lead to an underestimation of NUV SFRs 
while lower values lead to an overestimation in the NUV 
SFRs. Not only is Rv = 6 inconsistent with the value ob- 
tained for the matching between Ha and FUV SFRs, it is 
also too high to be a realistic value ( |Seaton|[l979| |Cardelli| 
et al.|1989||Calzetti||2001[). The answer to this dilemma lies 



in the | Cardelli et al.| (|1989j) curve. This obscuration curve, 
over-corrects the NUV derived SFRs more than the FUV 
derived SFRs. A similar result is obtained when the lSeatonl 



1 1979) obscuration curve is used. 



using the Cardelli et al. ( 1989 I curve for both the Ha and 



In comparing Figures [5] and [3] the Calzetti (2 001[ ) obscu- 
ration curve (used in panel (b) in both figures) is the only 
obscuration curve that does not over-correct NUV more than 
FUV. This is due to the inclusion of the 2200 A feature in 
the MW and FD05 extinction curves which magnifies the 
dust corrections. This feature is present in all the curves 
that are discussed here except for the Calzetti ( 2 00 1[ ) curve. 
The sudden increase in the extinction correction near NUV 
wavelengths is undoubtedly the source of the overcorrection 
seen in Figure [3] For this reason, we calculated the NUV 
dust correction using an FD05 curve with an R v value of 
4.5 but interpolating across the 2200 A feature (as shown 
in Table [T] where A v is the total attenuation). The results 
seen in Figure [3]d are very promising. The Ha and NUV 
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Figure 1. SFR derived from FUV luminosity as a function of SFR derived from Ha luminosity. All galaxies in our sample are shown, 
with AGNs excluded, (a) Both SFR indicators without obscuration corrections, and (b) only the Ha derived SFRs are corrected for dust 
obscuration using a FD05 curve with an R v value of 4.5. The axis ranges in the above panels, and subsequent figures, are the same for 
ease of comparison. The strong luminosity dependence of the obscuration corrections is clear (Hopkins ct al. 2001, 2003; Afonso et al. 
2003). 




SFR„„ (M.yr- 1 ) SFR Ha (M e yr _1 ) SFR Ha (M e yr" 1 ) 

Figure 2. FUV SFRs as a function of Ha SFRs, dust-corrected with a variety of obscuration corrections. Both Ha and FUV luminosities 
are corrected with consistent obscuration curves as described in the text, (a) A Cardclli et al. (1989) extinction curve is used to correct 
both Ha and FUV luminosities, (b) A Calzctti (2001) curve corrects both Ha and FUV luminosities, (c) Cardclli ct al. (1989) and 
Calzetti (2001) curves correct for the nebular and continuum emission respectively, (d), (e) and (f) show SFR values corrected using 
FD05 extinction curves with R v values of 3.5, 4.5 and 5.5 respectively for both luminosities. The best agreement between Ha and FUV 
derived SFRs is the FD05 extinction curve with an R v value of 4.5, panel (c). 
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to non-starburst galaxies the value of /3 will be overesti- 
mated, leading to an overestimation in the dust present in 
the galaxy and consequently the inferred SFR. 

Panel (b) of Figure 6 shows FUV derived SFRs cor- 
rected using a modified /3 attenuation correction to that of 
Meurer, Heckman and Calzetti (19991. The corrected at- 



tenuation factor was derived from the offset seen in FUV 
derived SFRs compared to Hen SFRs when the attenuation 
factor of Meurer, Heckman and Calzetti (1999) is used to 



correct FUV luminosities. We describe the attenuation as 

A = 3.3 4- 1.99/3. (15) 

The galaxies in Figure 6(b) do appear somewhat more evenly 
distributed on either side of the one-to-one line. There re- 
mains, though, a very large scatter, with an rms of 0.27. 

The conversion from luminosities to SFRs is carried out 
with the assumption that the IMF is constant and indepen- 
dent of the total SFR. Evolving IMFs could have an impact 
on the way luminosities are converted into SFRs, affecting 
the trends observed so far, but addressing this is deferred to 
later publications. 



SFRs are very closely matched, however, there appears to 
be a slight overestimation in the NUV SFRs at all SFRs. 
This is a result of the interpolation technique used to recre- 
ate this obscuration curve without the 2200 A feature. We 
find that to correct NUV luminosities for the effects of dust, 
the role of the 2200 A feature will have to be reassessed. The 
implication is that the 2200 A feature is a property of MW 
dust grains that does not seem to propagate to the global 
integrated attenuation properties of galaxies. 

Figure [4] shows a good correlation between Ha and [On] 
derive d SFRs. Both the FD05 (R v = 4.5) and |Cardelli et al.| 
( |1989[ ) curves show good agreement between Ha and [On] 
SFRs as seen in Figure [4] Figure [4]a indicates that [On] 
luminosities are obscured more than Ha as the Ha SFRs 
are larger than [On] before the dust corrections. This is ex- 
pected as shorter wavelengths are more sensitive to dust 
obscuration. We find that the metallicity dependence of the 
[On] SFRs does not play a significant role in biasing the 
SFRs derived from [On] for this sample although it may 
contribute to the observed scatter. It may become signifi- 
cant for the low-SFR systems which begin to deviate from 
the one-to-one line in Figure [4] ( K ewley et al.|2001 1 . 

Finally, we turn to an investigation of the utility of the /3 
parameter for making obscuration corrections to UV-derived 
SFRs. As mentioned above, /3 is demonstrated to be a useful 
obscuration parameter only for starburst galaxies (Kong et 
al. 2004). Our sample consists of galaxies with a broad range 
of SFRs, but containing very few starbursts. The distribu- 
tion of stellar-mass doubling-times, td, ranges from 10 s to 
10 13 years, with the vast majority (« 97%) having td > 10 9 . 
The definition of a starburst is such that their mass doubling 
times should be less than « lGyr, so even the very high SFR 
systems in our sample (in particular) are not starbursts, as 
they are also among the most massive. 

Our results (Figure 6) clearly show FUV SFRs, cor- 
rected for obscuration using f3, that are significantly overes- 
timated, by up to two orders of magnitude when compared 
to Ha derived SFRs corrected using the FD05 curve. As 
Kong et al. (2004) suggest, when the f3 parameter is applied 



5 SFR HISTORY 

The SFR distributions shown in the previous section include 
galaxies at various stages of their evolution. By comparing 
these galaxies with evolutionary synthesis models, using PE- 
GASE.2 ( |Fioc fc Rocca-Volmerange|p97t , we can develop 
an understanding of the effects that various galaxy param- 
eters have on their evolution and even explain some of the 
trends and scatter observed in our data. The output of the 
evolutionary synthesis model is normalized to a galaxy of 
one solar mass. To compare the theoretical models with the 
observed data we must scale the theoretical output accord- 
ing to the mass of the galaxies we are analysing. 

The galaxy mass is derived using the g-band and i-band 
colours and the absolute magnitude in the i-band (Mi). 



M _ 1 gl.04+0.9(g-i)-0.4M i 



(16) 



where g and i refer to the SDSS Petrosian magnitudes ob- 
served in the g and i bands ( Baldry et al.|20 04 Talor et al 
|2009| ). 

Figure [7] shows how the modeled evolutionary paths 
compare to the observed measurements. The paths in Fig- 
ure [t] have a |Baldry fc Glazebrook] ( |2003| > IMF and an ex- 
ponentially declining SFR with r = 90Myr. To give an in- 
dication of time on the evolutionary paths, the evolutionary 
paths begin at Myr and are displayed to « 1000 Myr in 
Figure [8] The evolutionary paths indicate that the FUV de- 
rived SFRs and Ha derived SFRs should strongly agree at 
all the SFRs that we have obtained. For the case of Ha and 
FUV SFRs this is consistent with our measurements. 

For the majority of the galaxies we compare, the masses 
of observed galaxies and those of the evolutionary paths 
seem to agree at the ages of 200 to 500 Myr. In Figure |SJ 
these ages on the evolutionary paths correspond, given the 
exponentially declining SF histories, to SFRs ranging from 
10 1 to 10 _2 Mq yr~ x . This is consistent across all mass scales 
in Figure [7| The masses of the modelled paths and the ob- 
served data do not agree at younger ages indicating that we 
do not observe these galaxies at the earliest stages of their 
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Figure 3. NUV SFRs as a function of Ha SFRs dust corrected with a variety of obscuration corrections. Both Ha and NUV luminosities 
are corrected with consistent obscuration curves, (a) A |CardeIli eTaLl|T989l l curve was used to correct both Ha and NUV luminosities, 
(b) A Calzetti (2001) curve corrects for both Ha and NUV luminosities, (c) Cardelli et al. (1989) and Calzetti (2001) curves correct 
for the nebular and continuum emission respectively, (d) A FD05 curve with an R v value of 4.5 and with the 2200 A feature is used to 
correct both Ha and NUV. (e) A FD05 curve with an R v value of 4.5 and without the 2200 A feature is used to correct both Ha and 
NUV. (f ) A FD05 curve with an R v value of 6 and with the 2200 A feature is used to correct both Ha and NUV. The best agreement 
between Ha and NUV derived SFRs is when the FD05 extinction curve with an R v value of 6, panel (f), and 4.5 without the 2200 A 
feature, panel (e), are used. This result is at odds with that of Figure (2) as the extinction curves are not consistent between FUV and 
NUV luminosities unless we remove the 2200 A feature. 
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Figure 4. [On] SFRs as a function of Ha SFRs. The [On] luminosities were corrected based on the extinction at the Ha wavelength, 
(a) compares the SFR indicators without any extinction corrections, (b) compares SFRs obscuration corrected with the Cardelli et al. 
(1989) obscuration curve, (c) compares SFRs obscuration corrected with the FD05 obscuration curve with an R v value of 4.5. 



Obscuration in galaxies and SFHs 9 



12 



10 



^ 8 

■< 



-- FD05 fl„=3.5 
-•- FD05fl r =4.5 
!;»■)«. 2.5 
" > Calzetti (2001) 
■ ■ Cardelli (1989) 




0.1 



0.2 



0.3 



0.4 

(pm) 



0.5 



0.6 



Figure 5. The Cardelli et al. (1989), Calzetti (2001) and FD05 
(R v = 3.5, 4.5, 5.5) obscuration curves. The 3 unbroken lines show 
the FD05 curves with the 2200 A feature removed. It must be 
noted that the Calzetti (2001) curve is the only curve intrinsically 
without the 2200 A feature. 



lives. Particularly the "hook" like path seen at the early 
stages (< 200 Myr) is not observed in our data and does not 
seem to constitute any of the scatter . At the older end, the 
masses of the models and the observed galaxies agree up to 
about 500 Myr. 

Figure [8]a shows that different IMFs can produce dif- 
ferent evolutionary paths particularly early in the life of a 
galaxy. The |Baldry fc Glazebrook| (|2003[> IMF begins at a 
higher SFR and using a |Baldry fc Glazebrook] p(jo"3| ) IMF 
will result in a galaxy at a particular age being observed to 
have a higher SFR compared to |Salpeter| ( |1955"| ) , Kcnnicutt 



( |1983[ ) or |Scalo| ( |1998[ ) IMFs. The modeled paths assume a 
galaxy mass of 10 a Mq with a 90 Myr characteristic decay 
time. All the IMFs share similar paths with the Baldry & 
|Glazebrook| ( |2003[ ) IMF showing a better fit to the observed 
trend than the other 3 IMFs. The tracks assuming Salpctcr 
( 1955[ ); |Kennicutt| ( |1983| and |Scalo|fl998[ ) IMFs follow sim- 
ilar trends at almost the same rate of evolution. The evolu- 
tionary paths for the IMFs match with the observed values 
at about the age of 130 Myrs. 

Also of interest is the lower SFR end where the evolu- 
tionary paths drift away from the one-to-one line and to- 
wards higher FUV-derived SFRs. This is clearly visible in 
the modeled evolutionary tracks. Even though our data does 
not extend below 10 _2 Mq yr~ l there is a mild but definite 
trend towards higher FUV SFRs compared to Ha SFRs ob- 
served at the low SFR end. By varying the characteristic 
decay times of the evolutionary paths, as in Figure [8]b, var- 
ious turn-off points can be modeled. A characteristic decay 




10" 3 10" 2 10 1 10° 10 1 10 2 10 3 

SFR Ha (M yr 1 ) 

Figure 6. The figure compares the Ha derived SFR, extinction 
corrected using a FD05 extinction curve (R v = 4.5), against FUV 
derived SFRs extinction corrected using the (5 parameter, (a) 
shows the FUV SFR corrected using the attenuation factor in 
Meurer, Heckman and Calzetti (1999). (b) is the same as (a) but 
with FUV SFRs corrected for dust using a modified attenuation 
factor to that of Meurer, Heckman and Calzetti (1999) chosen to 
fit the one-to-one line. There is a large scatter with an overesti- 
mation of the FUV SFRs. 



time of 90Myr seems to agree well with the observed values. 
Possible causes for this behaviour will be discussed in § 6. 

Figures [7] and [8] show good model agreement with the 
observed SFR distributions. It is clear that the evolution of 
galaxies with 10" 2 M Q yr' 1 < SFR < 10 2 M^yr' 1 can be 
modeled effectively. 



6 DISCUSSION 

6.1 Obscuration curves 

We have presented an analysis of the dust corrections in a 
sample of 31058 galaxies from the GAMA survey. We anal- 
ysed several techniques used to correct for the dust obscu- 
ration in galaxies. 

The Calzetti formalism using the Calzetti (2001) obscu- 
ration curve has been deemed to be one of the best for ap- 
plying obscuration corrections to FUV emission in starburst 
galaxies but we find that in our sample of non-starburst 
galaxies, it overestimates the UV SFR. We find that a theo- 
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Figure 7. FUV SFRs as a function of Ha SFRs compared with PEGASE evolutionary models of galaxies with various masses. The 
colour scale represents the measured masses of the galaxies. The regions where the modelled paths agree with the observed masses give 
an indication of the ages of the galaxies in each part of each curve. The evolutionary paths are for exponentially decaying SFRs with 
masses of 10 s MQ (red), 10 9 M© (black) and 10 10 M© (blue). All paths have a characteristic decay time of 90Myr and and a |Baldry fc| 
|GIazebrook| (2003) IMF. 




SFR„„ (MaJT- 1 ) SFR Hia (M . \i ) 



Figure 8. FUV SFRs as a function of Ho SFRs compared with PEGASE evolutionary models for different IMFs, panel (a), and charac- 
teristic decay times, panel (b). (a) |Salpeter| jl955) (blue), [Ke"nnicutt|(|1983) ( green), |ScaIo| (1998) (red) and |Baldry fc Glazebrook| ( |2003| > 
(black) IMFs are shown. All IMFs are scaled to a galaxy of mass 10 a M© and a characteristic decay time of 90 Myr. (b) Characteristic 
decay times of 30 Myr (black), 60 Myr (green), 90 Myr (blue) and 120 Myr (red). All paths were derived using a |Baldry &; G lazcbrook 
j2003) IMF scaled to a galaxy of mass 10 9 M0. 



retically modelled obscuration mechanism by FD05 provides 
SFRs well-matched between FUV and Ha. In addition, we 
use a Balmer decrement to derive the colour excess suffered 
by the stellar continuum, E(B — V) s tar, to complete our UV 
dust corrections. 

The effect of applying dust corrections is shown in Fig- 
ure [T] where the large increase in SFRs after applying dust 
corrections (to the Ha emission in this case) is a consequence 
of high star formation rate galaxies having much higher dust 
obscuration levels |Calzetti|2001| |Hopkins et al.|2001||2003 



Afonso et al. 20031 van den Bergh 20071. The same figure 



illustrates the lack of dust in low star forming systems as 
galaxies at the very low SFR end maintain similar SFRs 
after dust corrections. 

We explored several obscuration curves used to correct 
emission line luminosities, such as the |Seaton| ( |1979[ ) and 
Cardelli et al. ( 1989 1 MW obscuration curves as well the 



FD05 curves. It is important that the appropriate type of 
obscuration curve is applied for different types of emission. 
The FUV emission from galaxies contain contributions from 
older and lower mass stellar populations (Bell 2003) which 
would have gradually moved away from the clouds of gas 
that formed them, whereas Ha and [On] lines are emitted 
from the [Hn] regions ionised by hot young stars that are 
contained within those clouds of gas. For this reason, differ- 
ent absorption curves could potentially apply to stellar and 
nebular emission, highlighting the importance of using the 
correct obscuration curve for each type of emission. 

This is clearly illustrated in Figures [2]a and |3]a where 
using the |Cardelli et al.| (1989) curve for both Ha and UV 
overestimate the UV SFRs. The |Cardelli et al.|(l989) curve 
is appropriate for correcting nebular emission, and is more 
likely to over-correct for dust when applied to the UV stellar 
continuum. 



Obscuration in galaxies and SFHs 11 



The overestimation of the UV SFRs when the Calzctti 
( 2001 I curve is used (Figures [2]b and [2]c and |3]b and [3]c 
can have several explanations. As highlighted by Calzetti, 
(Calzetti 2001 1, this curve is only applicable to starburst 



galaxies. It was also based on a small sample of 37 star- 
burst galaxies using a universal extinction as opposed to a 
galaxy specific approach. It must be noted that a starburst 
galaxy is not the same as a high SFR galaxy. A starburst 
galaxy requires a short mass-doubling time compared to a 
non-starburst galaxy reflecting an atypically high SFR. Most 
of the galaxies in our sample have mass doubling times con- 
sistent with not being in starburst mode. For this reason, it 



is conceivable that the Calzetti (2001) obscuration correc- 



tions would overestimate the UV derived SFR of the galaxies 
in our sample. The processes that lead to starbursts such as 
interaction of galaxies, tidal shear in the solid body rota- 
tion and secular evolution of bars ( Larson fc Tinsley||1978 



Sanders et al.|1988||Kennicutt|1990l|Norman, Hasan fc Sell-| 
wood|| 1996| |Kenney, Carlstrom fc Young||1993[ ) would not 



necessarily occur in non-starburst galaxies and hence the 
spectral output can be different to a starburst galaxy. 

The FD05 extinction curve, with an R v value of 4.5 
is successful in providing matching Ha, FUV and [On] de- 
rived SFRs at all SFRs that we have tested (10~ 2 Mq yr' 1 < 
SFR < 10 2 Mq yr' 1 ) An R v value of 4.5 a grees with values 
prescribed by |Calzetti| ( [200T| . The FD05 curves are theoret- 
ical models based on the physical density of the turbulent 
structure and hence do not carry with them any of the con- 
straints that other curves have in their application. As the 
model deals only with the absorption and has nothing to 
do with the properties of the emission sources, the FD05 
curves can be applied to any given emission as they are not 
situation dependent. 

In the NUV, using the same FD05 curve as for the FUV 
(R v = 4.5) leads to an overestimation in the NUV SFRs 
compared to Ha derived SFRs. There is agreement when 
R v — 6 but this is inconsistent with other measured values, 
and with the value of R v found to reconcile the FUV, Ha 
and [On] measurements in our analysis. The fact that the 
NUV SFRs are overestimated more than the FUV SFRs in 
all MW obscuration curves and not in the Calzetti (2001) 
curve suggest that this effect is due to the 2200 A feature 
being included in the extinction curve. The NUV being the 
only SFR indicator that otherwise does not agree with the 
other indicators lends weight to this argument as none of 
those wavelengths fall within the range of the 2200 A feature. 
We show that the Ha and NUV derived SFRs agree very well 
if the 2200 A feature is eliminated in an FD05 curve with an 
R v value of 4.5. This provides us with a consistent extinction 
correction curve for all the SFR indicators. The role of the 
2200 A feature will require further analysis in future work. 

The use of the /3 parameter in estimating obscuration 
corrections for our sample of largely quiescent galaxies shows 
that intrinsic SFRs will be overestimated. Kong et al. (2004) 
attribute this to a secondary effect based on SFR histories. 



6.2 SFR Histories 

We applied galaxy evolutionary synthesis models to the com- 
parison between Ha and FUV derived SFRs, and there was 
good agreement for all SFRs. The evolutionary paths were 



compared according to galaxy mass, characteristic decay 
time and IMF (assuming an exponentially decaying SFR). 

We showed that by comparing the evolutionary paths 
for varying galaxy masses with observed data, we can iden- 
tify a range of ages that the observed galaxies would fall 
into. This range was estimated to be between 200 Myr and 
500 Myr. There is a continuous variation in the ages mak- 
ing it difficult to identify any effects on the ages of galaxies 
due to mass. The IMFs we have tested show similar trends 
with a |Baldry fc Glazebrook| ( |2003[ ) IMF showing the best 
agreement with the observed data (Figure [8|a). This is par- 
ticularly clear at the low SFR end with the |Baldry fc Glaze-| 
|brook| {2003} IMF following t he observed SFR trend better 
than the other 3 IMFs. The |Baldry fc Glazebrook| (|2003 \ 



IMF has a faster evolution, where at any given SFR, the 
|Baldry fc Glazebrook| ( |2003[ ) evolutionary path has an age 
that is about 50Myr higher than the other tested IMFs. 

The characteristic decay time was used to show how 
galaxies behave as they age: At low SFRs, a turn-off is ex- 
pected away from the one-to-one line, towards the FUV SFR 
axis, as the Ha derived SFRs decrease faster than FUV de- 
rived SFRs. Our data does not extend to such low SFRs, but 
using PEGASE models we can recreate a predicted turn-off 
had our data extended to such low SFRs. The turn-off point 
is highly dependent on the characteristic decay time of SFRs 
with higher characterised decay times causing the turn off 
to occur later. 

One of the proposed mechanisms for this behaviour 



is the Integrated Galaxy IMF (IGIMF) theory ( Pflamm- 
lAltenburg, Weidner & Kroupal 120071 120091). Pflamm- 



Altenburg, Weidner & Kroupa ( 2007 2009 1 constructed an 



integrated galactic initial mass function by combining all 
IMFs of all young star clusters. Using the IGIMF they pre- 
dicted that the Ha luminosity of galaxies should decrease 
faster than UV luminosities at low SFRs as the Ha lumi- 
nosity of star-forming dwarf galaxies decreases faster with 
decreasing SFR than the UV luminosity at low SFRs. This 
is clearly observed in the evolutionary paths in Figures |8]a 
and b. |Pnamm-Altenburg, Weidner fc Kroupa| p007] |2009 1 



predict this turn-off to occur at around 10~ Mq yr~ . As 
seen in Figure[8]b this turn-off can be replicated by a 90 Myr 
characteristic decay time which also agrees well with the ob- 
served data at higher SFRs. Even though our data does not 
extend far below an SFR of 10 _2 Mq yr -1 , the beginnings 
of such a trend can be observed in Figures [2]e,[7]and|8] 



7 CONCLUSION 

We find that an FD05 obscuration curve with an R v value 
of 4.5 gives the best agreement across all SFR indicators 
provided that the 2200 Afeature is not included. These SFR 
distributions agree well with modeled evolutionary paths. 
We also show that the /3 parameter is not suitable for use 
with GALEX data. 

Our main objective has been to obtain accurate dust ex- 
tinction corrections for Ha and FUV and present a clear and 
concise formalism on how this is to be done based on various 
obscuration laws described in the literature. While the dust 
corrections were largely successful, we have identified key 
areas for further research, particularly in the significance of 
the 2200 A in obscuration corrections. Analysing the effec- 
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tiveness of extinction curves, particularly with IR data will 
be of great significance for future work. 
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